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Abstract. - To clarify the order parameter symmetry of cuprates, the magnetic penetration depth 
A was measured along the crystallographic directions a, b, and c in single crystals of YBa2Cu408 
via muon spin rotation. This method is direct, bulk sensitive, and unambiguous. The temperature 
dependences of A^" 2 and A^ 2 exhibit an inflection point at low temperatures as is typical for 
two-gap superconductivity (TGS) with s + d— wave character in the planes. Perpendicular to the 
planes a pure s-wave gap is observed thereby highlighting the important role of c-axis effects. We 
conclude that these are generic and universal features in the bulk of cuprates. 



Two gap superconductivity (TGS) remained a theoret- 
ical issue only for more than 20 years [1—3], even though 
it seemed to be a natural and intriguing extension of BCS 
theory for more complex materials. In 1980 TGS was fi- 
nally observed in Nb doped SrTi03 [4] and believed to be a 
rare exception in superconductors. Since the discovery of 
TGS in MgB2 many more superconductors with TGS were 
found, including heavy fermion compounds, making this 
feature more common than believed early on. In all above 
mentioned superconductors the combined order parame- 
ters exhibit always the same symmetry, namely s + s in Nb 
doped SrTiOa [4] and MgB2 [5] and d + d in heavy fermion 
compounds [6] . A totally novel situation is met in cuprate 
high-temperature superconductors (HTS), since the order 
parameters are of different symmetries, i.e., s + d [7-13]. 
However, theoretical modelling suggested a single d— wave 
order parameter in the Cu02 planes and unfortunately 
biased further research and partly inhibited the experi- 
mental efforts to characterize TGS in more detail in HTS. 

Recently, new muon spin rotation (/iSR) investiga- 
tions of single crystal Lai. 8 3Sr .i7CuO4 detected an inflec- 
tion point in the temperature dependence of the inverse- 
squared in-plane magnetic penetration depth A~ fc 2 , which 
is a direct consequence of two superconducting gaps with 



largely different zero temperature gap values [14]. Since it 
is unclear whether these observations are a material spe- 
cific property or generic and intrinsic to HTS, the previous 
/iSR studies were extended to single crystal YBa2Cu40s 
and were performed by applying a magnetic field along the 
crystallographic directions a, b, and c. Thereby we obtain 
the three principle components of the second moments of 
the local magnetic field distribution P(B) in the mixed 
state, which are related to the superfluid density, and re- 
flect directly the corresponding penetration depths A a , A;,, 
and A c . While A~ 2 and \Z 2 vary almost linearly with tem- 
perature for 20 < T < 50 K, as is expected for a d— wave 
order parameter, the c— axis response (A~ 2 ) saturates be- 
low 30 K, as expected for a s— wave order parameter. In 
addition, A~ 2 and A^" 2 , both exhibit an inflection point 
in their temperature dependences around T ~10 K which 
- as has been shown in [15] - is the consequence of two 
coexisting order parameters, namely s + d. 

Details of the sample preparation for YBa2Cu40s can 
be found elsewhere [16]. All crystals used in the present 
study were taken from one batch. The superconducting 
transition temperature T c and the width of the supercon- 
ducting transition AT C were determined for the three sets 
(~40-50 crystals each) of the main set (~130) of crys- 
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tals. Both were obtained from field-cooled (0.5 mT) mag- 
netization curves measured by a SQUID magnetometer 
and exhibited all the same values, i.e., T c ~ 79.9 K and 
AT C ~ 2 K. The crystals had mostly a rectangular shape 
with a typical size of approximately 0.8x0.3x0.05 mm 3 . 
X-ray measurements revealed that the crystallographic 
6— axis is exactly parallel to the longest side. Bearing in 
mind that the c— axis is perpendicular to the flat surface 
of the crystal, we were able to orient the whole set along 
the crystallographic a, b, and c directions. The final orien- 
tation of the crystals in the mosaic was checked by using 
a polarizing microscope. 

The transverse-field /iSR experiments on a mosaic of ori- 
ented YBa2Cu408 single crystals were done at the 7rM3 
and 7rEl beam lines at the Paul Scherrer Institute (Vil- 
ligcn, Switzerland). The mosaic was field cooled from 
above T c to 1.7 K in a field of 0.015 T. The /xSR ex- 
periments were performed for the magnetic field applied 
parallel to the a, b, and c crystallographic axes. Typical 
counting statistics were ~24-30 million muon detections 
over three detectors. In the analysis presented below we 
used the well-known fact that for an extreme type-II su- 
perconductor in the mixed state A -4 is proportional to 
the second moment of P{B) probed by /xSR [17]. The 
second moment of P{B) was calculated within the same 
framework described in [14]. We used a four component 
Gaussian expression to fit the /iSR time spectra. One com- 
ponent arises here from the background signal stemming 
from muons stopped outside the sample, whereas the other 
three components describe the asymmetric line shape of 
P(B) in the mixed state (see inset to Fig. 1). The first and 
the second moments of P(B) (excluding the background 
component) are then obtained as: 
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Here 7^ = 2ir x 135.5342 MHz/T is the muon gyromag- 
netic ratio. Ai, (Tj, and Bi are the asymmetry, the relax- 
ation rate, and the first moment of the i— th component, 
respectively. The superconducting part of the square root 
of the second moment (<r sc oc A~ 2 ) was then obtained by 
subtracting the nuclear moment contribution (a nm ) mea- 
sured at T > T c , according to <r 2 c = a 2 — er 2 m . To ensure 
that the increase of the second moment of the measured 
/xSR signal below T c is attributed entirely to the vortex 
lattice, zero-field /iSR experiments were performed. No 
evidence for static magnetism in the YBa2Cu40s mosaic 
sample down to 1.7 K was observed. 

For an anisotropic London superconductor the effective 
penetration depth for the magnetic field along the i-th 
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Fig. 1: (Color online) Temperature dependences of a a i, oc A ai) 2 
(H ||c), a ac oc A~ c 2 (H\\b), and Ub c oc 

Kc ( H \\ a ) ofYBa 2 Cu 4 8 , 
measured after field cooling the sample in fioH =0.015 T. The 
inset shows the local magnetic field distribution P(B) obtained 
by means of the maximum entropy Fourier transform technique 
at T = 1.7 K and fi^H = 0.015 T applied parallel to the c— axis 
(open circles: data; solid line: four component Gaussian fit 
excluding the background). 



crystallographic axis is given by [18]: 
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Here the index "sc" for the superconducting part of the 
square root of the second moment a sc is omitted for sim- 
plicity. From Eq. (3) it is obvious that for the magnetic 
field applied along one of the principal axes a, 6, and c, 
the components Obc oc A^ 2 , <j ac oc A~ c 2 , and a a b oc A~ 2 
are measured. The temperature dependences of <7b c , <7 ac , 
and <j a b after field-cooling the sample in fj, H = 0.015 T 
are shown in Fig. 1. It is seen that at Ti p ~ 10 — 20 K 
all measured <Jij(T) exhibit an inflection point. Below 
this point a a b, <7bc, and a ac increase by approximately 
70%, 35%, and 10%, respectively. Note, that a similar 
inflection point in A~ h 2 (T) was also observed in low-field 
magnetization (LFM) experiments on powder samples of 
YBa2Cu408 [19]. However, in these experiments the in- 
crease of A~ 6 2 below Tip was much less pronounced than 
the one observed in the present study. This difference can 
be explained by the fact that LFM probes the penetration 
depth mainly near the surface, whereas fiSH measures A 
in the bulk. In LFM experiments the magnetic field pene- 
trates the sample only on a distance A from the surface of 
the sample (few hundred nanometers), thereby leaving the 
main part of the superconducting volume unaffected. In 
contrast, muons penetrate at much longer distances into 
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the sample (few-tenth of a millimeter), thus probing A 
deeply inside the sample. It should be recalled that at 
the surface due to the symmetry breaking, the order pa- 
rameter has pure d— wave character, whereas in the bulk 
s + d— wave symmetry is present [11]. This explains in a 
natural way the difference between the surface sensitive 
LFM and the bulk sensitive /iSR experiments. 

From Eq. (3) the individual components of A along the 
i— th crystallographic direction are obtained as: 



Figure 2 shows the temperature dependences of a a cx A~ , 
<7(, cx A^ 2 , and a c cx A~ 2 as derived from the data pre- 
sented in Fig. 1 by using Eq. (4). In the following we dis- 
cuss the temperature dependence of each particular com- 
ponent separately. 

From Figs. 2 (a) and (b) it is evident that both <r a cx A~ 2 
and <7b cx A^ 2 increase almost linearly with decreas- 
ing temperature in the range 50 < T < 20 K. Around 
T ip ~ 10 K an inflection point is visible in the temper- 
ature dependence. Below this point an unusual increase 
in both quantities appears: <j a increases by almost 20% 
(from 2.95 ^s" 1 at T = 10 K to 3.5 ^s" 1 at T = 2.6 K) 
and a b by 60% (from 5.2 ^s" 1 at T = 10 K to 8.3 /is" 1 at 
T = 2.6 K). As has been shown in Ref. [15], an inflection 
point in A~ 2 (T) suggests the presence of at least two su- 
perconducting gaps in YBa2Cu40g with very different gap 
values, i.e., a large gap and a small one. The same behav- 
ior was observed recently in Lai.83Sro.i7Cu04 by /iSR [14], 
as well as in other HTS using various techniques [9-13], 
supporting a two-gap behavior with the larger gap being 
of d— wave and the smaller one of s— wave symmetry. The 
above data were analyzed by assuming that an s— wave 
and a d— wave gap contribute to a according to: 

a(T) = a d (T)+<7 s (T), (5) 

where both components are expressed like [14]: 

g(T,A ) =1+ 1 ^ r (df_\ EdEdy 

(6) 

Here, / = [1 + exp(_E/fcsT)] 1 is the Fermi function, 
Ao is the maximum value of the gap, and A(T,tp) = 
A A(T/T c )g((p). For the normalized gap A(T/T C ) tabu- 
lated values of Ref. [21] were used. The function g(ip) de- 
scribes the angular dependence of the gap and is given by 
g d (ip) = | cos(2y>)| for the d— wave gap [12] and g s (tp) = 1 
for the s— wave gap [see insets in Figs. 2 (a) and (b)]. 
From this analysis we obtain Ag a = 19.41(9) meV, <r~ = 
3.03(2) us- 1 , A s Q a = 0.96(5) meV, a s a = 0.48(2) ^s" 1 
and A d b = 21.17(8) meV, a d = 4.66(2) ^s^ 1 , Ag b = 
1.10(2) meV, a s b = 3.63(2) ^s" 1 for a a (T) and a b {T), re- 
spectively. The main results are summarized in Table 1. 

The temperature dependence of er c cx A,T 2 differs sig- 
nificantly from the one of the in-plane components since 
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Fig. 2: (Color online) Temperature dependences of a a cx A„ 2 
(a), at cx (b), and a c cx A~ 2 (c) of YBa2Cu408 at 

fioH =0.015 T obtained from the data presented in Fig. 1 
by means of Eq. (4). Lines in (a) and (b) represent results 
of the fits of <Ja{T) and o"b(T) by means of the two-component 
[Eq. (5)] and the three-component models. The blue dotted 
curve in (c) represents a fit to the data assuming that cr&(T) is 
entirely determined by the in-plane d— wave and the chain gap 
measured in [20] (see text for details) . The black solid line and 
the blue dotted curve in (c) represent fits by means of Eq. (6) 
with g 3 {ip) = 1 and the power law a c (T) = cr c (0)[l - (T/T c )"], 
respectively. The insets in (a) and (b) represent the angular 
dependences of the d— wave gap (A d ), the s— wave gap (A s ) 
and the gap in CuO chains (A ch ) obtained in [20]. The inset 
in (c) shows a c (T) at T < 30K. 

it saturates at temperatures T < 30AT to become T in- 
dependent [see inset in Fig. 2 (c)]. This dependence is 
analyzed in two ways: (i) by using Eq. (6) with g s (<p) = 1, 
yielding Ag c — 19.20(4) meV, and (ii) by assuming the 
phcnomenological power law u c (T) = <r c (0)[l — (T/T c ) n ]. 
The results are given in Table 1 and compared to the ex- 
perimental data in Fig. 2 (c), from which it is obvious 
that both approaches are almost undistinguishable. For 
the power law dependence a critical exponent n — 4.5(3) 
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Table 1: Summary of the gap analysis of the temperature de- 
pendences of a; oc A~ 2 (i — a, b, c) for YBa2Cu40s. The 
meaning of the parameters is explained in the text. 



Component 


\d 

(meV) 


(meV) 


lab 


lab 


<Ja 


19.41(9) 


0.96(5) 






Ob 


21.17(8) 


1.10(2) 


1.09(1) 


1.15(6) 


Oc 




19.20(4) 







is derived which is close to the one obtained within a two- 
fluid model where n = 4, which applies to a strong cou- 
pling s— wave BCS superconductors [22]. Using a d— wave 
model for a clean superconductor, Ref. [23] predicts an 
exponent n = 5. Since the observed power law expo- 
nent lies between these two values, it could be argued 
that d— wave superconductivity could also be the cause 
of the temperature dependence of a c oc A~ 2 . However, 
this can be ruled out, since the results of Ref. [23] do 
not apply to cuprate superconductors containing chains, 
like YBa2Cu307_a and YBa2Cu40s- We can thus safely 
conclude that the order parameter along the c— axis is 
isotropic s— wave. This conclusion is further supported 
by c— axis tunneling [24], bicristal twist Josephson junc- 
tions [25] , optical pulsed probe [26] , and optical reflectivity 
(see Fig. 4 in Ref. [27]) experiments, as well as theoretical 
considerations [28] . It is important to note that due to the 
reasons described above the saturation of A~ 2 (T) obtained 
here has not been observed in LFM experiments [19]. 

The anisotropy observed in a a and <7b and their un- 
conventional temperature dependences deserve further re- 
marks. The temperature dependences of both quantities 
are characterized by an inflection point at Ti p ~ 10 K evi- 
dencing that at least two superconducting gaps contribute 
to the superfluid density. However, along the crystallo- 
graphic b— direction this is much more pronounced than 
along a, and the analysis of both data sets in terms of 
s + d wave gaps using Eq. (5) reveals that along the a 
direction the s— wave order parameter contributes 14%, 
whereas along the b direction the s— wave contribution is 
42%. This observation is a consequence of the structural 
anisotropy caused by the CuO chains constitutes itself by 
an enhanced s— wave superconducting density along the 
chain direction which can have various origins. Early on, 
Raman experiments have detected the opening of a second 
gap at low frequencies which appears only in specific scat- 
tering geometries [29]. Since these observations have only 
been made in the YBCO family, they were interpreted as 
being due to the opening of a gap in the chains. In addi- 
tion, recent angle resolved photoemission data (ARPES) 
have observed a coherence peak related to the chain bands, 
which appears in a very limited k— space region [20]. It 
was found that in YBa 2 Cu 4 08 the superconducting gap 
in the chains opens only at certain angles if and reaches 
its maximum value of 5(1) meV at <p ~ 60° [see inset in 



Fig. 2 (b)]. Therefore, one may speculate that the two- 
gap behavior in YBa2Cu 4 08 is entirely determined by the 
in-plane d— wave and the chain gap. This scenario can, 
however, be excluded. First of all, the fit of Eq. (5) to 
Ob (T) taking into account a chain gap of this symmetry, 
leads to a rather poor agreement with the experimental 
data [blue dotted curve in Fig. 2 (b)]. Second, the ob- 
servation of an inflection point in cr a (T) cannot be ex- 
plained by a possible misalignment of the crystals in the 
mosaic. <7*(0) = 0.48 [is^ 1 would correspond to a sit- 
uation when all the crystals in the mosaic are turned by 
±8°, or ~30° FWHM if one assumes a triangular distribu- 
tion of orientations with the maximum along the crystallo- 
graphic b axis. Therefore, the data presented in Fig. 2 (b) 
were analyzed by assuming that &b(T) is the sum of three 
components: cr^(T), <rg(T), and the contribution from the 
chains a^ h (T). These results are, however, undistinguish- 
able from the s + d analysis presented above. It is impor- 
tant to note, that independent of a possible chain related 
energy gap, both procedures provide clear evidence that 
the in-plane superconducting order parameter consists of 
two components, a d— wave order parameter and an addi- 
tional s— wave component. 

The maximum values of the in-plane d— wave gap along 
the a and b directions are of similar order of magnitudes, 
i.e., A^ a = 19.41(9) meV, A£ 6 = 21.17(8) meV and thus 
in a good agreement with Aq = 22 meV derived from tun- 
nelling experiment [30]. The in-plane anisotropy in both 
gap values, s and d, is approximately the same, i.e., 7^ = 
Ag i6 /Ag >0 - 1.09(1) and 7a A ; - A^/Ag^ = 1.15(6) (see 
Table 1). Moreover, it is important to note that the 
s— wave gap along the c direction Aq c — 19.20(4) meV 
is of the same order of magnitude as the d— wave gaps in 
the ab plane: Ag c ~ Aq a ~ A£ 6 « 20 meV. 

In conclusion, we performed a systematic /xSR study 
of the magnetic penetration depths A a , Ab, and A c on 
single crystals of YBa2Cu40s. In contrast to previous 
LFM experiments [19], our method is bulk sensitive and 
a direct probe of the penetration depth. The use of sin- 
gle crystals enables us to derive the magnetic penetra- 
tion depth along the three principal crystallographic di- 
rections. Along the a and b directions clear evidence is 
obtained that TGS is realized also in the chain containing 
compound YBa 2 Cu40s. While the in-plane penetration 
depth is anisotropic, exhibiting an inflection point at low 
temperatures in both A~ 2 and A^" 2 which is characteris- 
tic of TGS, the c— axis data provide clear evidence for 
an an isotropic s— wave gap. From the data it must be 
concluded that the in-plane superconducting gap consists 
of two components, namely s + d. The situation along 
the c— axis is different and supports s— wave only. This 
exceptional behavior has not been predicted by any the- 
ory, since the third dimension has mostly been neglected. 
Since YBa 2 Cu40s differs structurally substantially from 
the previously investigated system L^-zSr^CuO^, where 
the in-plane penetration depth also shows s+d wave super- 
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conductivity [14] , the above findings cannot be attributed 
to specific structural features of the chain containing com- 
pound, but show that s + d order parameters are generic, 
intrinsic and common to all HTS. The results thus exclude 
theoretical approaches as e.g., the "plain vanilla" mecha- 
nism which concentrates on the Cu0 2 plane only. Besides 
of the fact that t-J or 2D Hubbard models are not capa- 
ble to yield the observed s— wave component, the role of 
the lattice played for superconductivity has to be reconsid- 
ered since the presence of an s— wave component naturally 
points to its importance. 
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